Block scaling in the directed percolation 
universality class 



Gunnar Pruessnerf X 

f Mathematics Institute, University of Warwick, Gibbet Hill Road, Coventry CV4 
7AL, United Kingdom 

^Present address: Department of Mathematics, Imperial College London, 180 
Queen's Gate, London SW7 2AZ, United Kingdom 

E-mail: g . pruessner@imperial .ac.uk 



Submitted to: NJP 

PACS numbers: 05.70.Jk, 05.40.-a, 89.75.-k 

Abstract. The universal behaviour of the directed percolation universality class is 
well understood, both the critical scaling as well as finite size scaling. This article 
focuses on the block (finite size) scaling of the order parameter and its fluctuations, 
considering (sub-)blocks of linear size / in systems of linear size L. The scaling 
depends on the choice of the ensemble, as only the conditional ensemble produces 
the block-scaling behaviour as established in equilibrium critical phenomena. The 
dependence on the ensemble can be understood by an additional symmetry present in 
the unconditional ensemble. The unconventional scaling found in the unconditional 
ensemble is a reminder of the possibility that scaling functions themselves have a 
power-law dependence on their arguments. 
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1. Introduction 

The directed percolation (DP) universality class comprises a huge number of non- 
equilibrium critical phenomena. Janssen and Grassberger jH [2] famously conjectured 
more than 25 years ago that under very general circumstances, all models with a unique 
absorbing state belong to the DP universality class. While this conjecture has been 
confirmed numerically many times, evidence for the presence of the DP universality 
class in natural systems is still very scarce [3] (but see [1]). 

One problem, when probing field data for the presence of DP is that field data is 
more readily obtained in a single measurement rather than as a time series. However, 
the statistical features to be identified require an entire ensemble of realisations of the 
process in question. Instead of using a time series, one can resort to sub-sampling, i.e. 
splitting a large sample of size L d into {L/l) d small ones of size l d . For example, a 
population pattern obtained by measuring the spatial distribution of species could be 
split into several distinct blocks and their mutual correlations analysed. 

The question how the order parameter of an equilibrium system at the critical 
point changes with the block size it is averaged over, has been studied in great detail by 
Binder [5]. In the present work, a corresponding analysis is applied to models belonging 
to the DP universality class, more specifically to the contact process and to directed 
percolation itself from the point of view of absorbing state (AS) phase transitions. It 
turns out that the block averaged order parameter needs to be defined very carefully in 
order to reproduce standard finite size scaling. In the following, I will present numerical 
evidence and theoretical arguments that block finite size scaling (FSS) in DP can be 
very different from what is expected from equilibrium critical phenomena, depending 
on the choice of the ensemble. This observation can be readily applied to the analysis 
of field data, and will be illustrated using surrogate data. 

2. Method 

The order parameter of an absorbing state phase transition is the activity p, which, 
in lattice models, is the density of active sites. The activity vanishes as soon as the 
system hits the absorbing state, from where it cannot escape. In a finite system, the 
absorbing state is reached with finite probability from anywhere in phase space, so that 
every (finite) system eventually becomes inactive, lim^oo p(t) = 0, where t measures 
the time in the model. However, the order parameter p signals a phase transition 
in a temperature-like tuning parameter p, dividing the parameter space into a region 
where the decay of the activity p is exponentially fast in time, from a region where for 
sufficiently large systems it is practically impossible to observe p = 0. 

At least two (seemingly) different methods have been devised to overcome the 
problem that strictly lim^oo p(t) = and to obtain the phase transition even in finite 
systems: One either introduces an external field which induces activity in the system [6] 
and analyses the model in the limit of arbitrarily small fields, or one considers the quasi- 
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stationary state [7J. In this latter approach, after initialising and discarding a transient, 
all averages are taken conditional to activity. The activity entering the observables 
therefore never vanishes and the order parameter is always non-zero, rendering for 
example moment ratios well defined. One can show that both methods produce 
asymptotically equivalent scaling results [8]. In the present article, the quasi-stationary 
state was used to produce individual samples. 

Above the critical point p c , i.e. for p > p c , the ensemble average of the activity p 
in the thermodynamic limit picks up as a power-law (p) = A(p — Pc) 13 , where A is the 
amplitude, p c is the critical value of the tuning parameter and (3 is a universal critical 
index. Similarly, fluctuations of the order parameter scale as a 2 (p) = (p 2 ) — (p) 2 = 
B±L~ d \p — p c |~ 7 , where B + and B_ are the amplitudes above and below the transition 
respectively, L is the linear extent of the system, d is its spatial dimension and 7 is an 
independent critical exponent. Both these power laws are asymptotes and thus acquire 
corrections [9] away from the critical point. 

Ordinary FSS pHl EED, EE] is observed by tuning p = p c and considering the 
behaviour of the observables as a function of the system size, (p) (L) = A'L" 13 ^ 1 - and 
cr 2 (p) (L) = B'L~ d+1 l Ux - , where v±_ is a third critical exponent. These exponents are 
related by — d + r y/v± = —2j3/b , ± [13J, which in equilibrium corresponds to Josephson 
hyper-scaling together with the Rushbrooke scaling law. 

Another kind of FSS can be explored in addition to the ordinary one just described. 
Instead of considering the entire system, observables are recorded within small blocks 
of linear extent /. To this end, I introduce the observables (p) u (I; L) and a 2 (p) (/; L), 
which are first and second cumulants of the activity within those blocks of linear size 
I in a system of size L. The blocks are produced by dividing samples of linear size L 
generated in the quasi- stationary state into (L/l) d blocks of linear size I each. This 
length I is conveniently chosen so that it divides L. By construction, each sample of 
size L d contains at least one active site and consequently at least one block of size l d 
has non- vanishing activity. However, there may be up to {L/l) d — 1 inactive blocks. 

In addition to these observables, I introduce (p) c (/; L) and a 2 (p) (I; L), which are 
the first and second cumulant of the activity conditional to activity within the respective 
block. That means that inactive blocks are discarded when averaging. To distinguish 
the two ensembles, the former (with cumulants (p) u (I; L) and a 2 (p) (/; L)) will be called 
"unconditional" (subscript u) and the latter (with cumulants (p) c {I; L) and a 2 (p) {l;L)) 
"conditional" (subscript c). One can derive the moments of the unconditional ensemble 
from the corresponding moments in the conditional ensemble and vice versa by re- 
weighting, because they differ only by a number of samples with vanishing block-activity 
which are discarded in the conditional ensemble but not in the unconditional one (see 
[8] for similar considerations for the overall activity). If the fraction of active blocks 
averaged over the entire unconditional ensemble is a(l; L), then the nth moment of the 
unconditional ensemble (p n ) u is related to the nth moment of the conditional ensemble 



(p n ) c by 



( P n ) u (I; L) = a(l; L) (p") c (/; L) 



(1) 
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for n > 0. 

More than 25 years ago, block scaling was investigated by Binder j5] for the Ising 
model, where the order parameter is the magnetisation density m rather than the 
activity p. In these systems with a symmetric phase-space, there is no corresponding 
distinction of active and inactive blocks. Naively transferring these results to DP 
suggests 

(p) u>c (l;L) = C u J-^g u , c (l/L) (2) 

in leading order of L and in the limit of / being large compared to a lower cutoff Iq, 
i.e. I ^> Iq. Below this constant threshold Iq, i.e. I < l , the order parameter deviates 
from the behaviour predicted by (j2J). This phenomenon is known and well understood 
in classical critical phenomena [T4"l [12] . In the following, data for I < l is not shown; 
in the two-dimensional contact process and directed percolation l was estimated to be 
Iq ~ 8 and in the one-dimensional contact process Iq ~ 16. 

The dimensionless scaling functions G u ,c( x ) i n (ED are bounded from above and 
usually also from below (away from 0) for all x G [0, 1]. The scaling function is universal 
up to a pre-factor, which can always be absorbed into the metric factor C UjC . The latter 
is required for dimensional consistency. 

Scaling behaviour similar to §Z§ is expected to hold for higher order moments and 
cumulants as well, in particular for the variance of the activity 

< c (p) (I; L) = D u J- 2 ^^ c (l/L) , (3) 

using the scaling relations cited above. Again, the scaling functions T UjC are constraint 
by being bounded from above and (usually) away from zero from below. 

So far, I have described what is expected in absorbing state systems as inferred from 
equilibrium critical phenomena. In the remainder of this article, I will first introduce 
the non-equilibrium models and the methods used in the numerical simulations for this 
study. The results for the different observables (activity, its variance and a moment 
ratio) in the different ensembles (unconditional and conditional) are then discussed in 
the light of analytical arguments. The article finishes with an application to surrogate 
data, a discussion of the implications and the wider context of the findings and concludes 
with a brief summary. 

3. Results 

Most of the numerics in this work is based on the contact process [TBI [T3] on a two- 
dimensional square lattice with periodic boundary conditions, but the same results are 
found for the one-dimensional contact process as well as for two-dimensional site-directed 
percolation [16J. In fact, it will be argued that they are general features of the directed 
percolation universality class, if not of all absorbing state phase transition. In the two- 
dimensional contact process occupied sites turn empty with extinction rate 1, while 
empty sites become occupied with rate zX with z being the fraction of occupied nearest 
neighbours. The critical value of A has been estimated numerically with great accuracy 




Figure 1. Data collapse of the order parameter measured in the unconditional (open 
symbols) and the conditional ensemble (filled symbols) for the two-dimensional contact 
process with system sizes L = 32, 64, . . . , 512 (up to L = 1024 in the conditional 
ensemble) and I = 16, 32, . . . , L (data points for the same system size obtained 
from the same samples are connected by lines). The unconditional order parameter 
collapses trivially (see text) and is shown for illustration purposes only. The sloped 
dashed line has slope 0.795 « indicating the power-law behaviour of the scaling 

function of the unconditional ensemble, Q u , defined in ([2]). The conditional order 
parameter collapses well under the scaling ansatz ([2]), plotting (p) c (I; L)l^/ U± vs. l/L. 
The horizontal dashed line is the likely asymptote of the scaling function Q c in the 
conditional ensemble. 



A c = 1.64877(3) jTTJ (I used A = 1.6488 in the present study). The time scale is set by 
the extinction rate. The two-dimensional contact process belongs to the DP universality 
class which is characterised by exponents (3 = 0.583(3), 7 = 0.297(2) and v± = 0.733(4) 
PEE], so that /3/v± = 0.795(6). 

In site-directed percolation in 2 + 1 dimensions (BCC lattice) the time evolves 
discretely [T6l [19] . A site is occupied in the following time step with probability p if 
at least one of its directed neighbours is occupied, otherwise it is empty The directed 
neighbours of a site are four sites in the preceding time step: The site itself, its right 
and upper nearest neighbour and its upper right next nearest neighbour. The critical 
value of p in this model has been estimated as p c = 0.34457(1) [18]. This model belongs 
to the 2D DP universality class as well. 
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All measurements are taken at the quasi- stationary state: Starting from random 
initial configurations with a small but non-vanishing activity, the systems evolve 
according to the rules described above, until the observables reach a (quasi-)stationary 
state. For example, in the two-dimensional CP for L = 256, the first 10 5 updates are 
discarded as transient. Measurements are taken at constant rate after the transient and 
enter with the same weight until the system reaches the absorbing state (for L = 256 the 
average lifetime was about 4.8 • 10 4 ) or a maximum time is reached (5 ■ 10 5 for L = 256). 
The procedure is repeated until the statistical error is acceptably small; for example 
1.42 • 10 s systems of size L = 256 were started, of which only about 13% survive the 
transient. Statistical errors were estimated by sub-sampling the ensemble, which copes 
even with correlated data. In the following, the scaling of the various (block) observables 
in I and L is analysed. 

3.1. Order parameter 

The first moment of the activity in the unconditional ensemble, (p) u (I; L), does not vary 
in I at all, because of translational invariance: Every (randomly chosen) site is equally 
likely to be active and therefore (p) u (I; L) is constant in I. In order to establish ordinary 
FSS complying to ([2]), the scaling function Q u {x) must necessarily be a power- law itself, 
Qu{x) = G u (l)x^ u , so that 

(p)Jl;L) = C u g u (l)L-^ (4) 

otherwise standard finite size scaling, (p) u (I = L; L) oc L~^l v ', would not be recovered 
for I = L, i.e. when a single block covers the entire system. Contrary to what is expected 
from equilibrium critical phenomena, this scaling function necessarily vanishes at 0, i.e. 
lim^o Gu(x) = 0. Figure [H contains a (trivial) data collapse for (p) u (/; L) according to 
(EJ), namely (p) u (l;L)0^ oc (l/L)^ function oil/L for various system sizes L, 
which illustrates the scaling behaviour. 

The conditional order parameter (p) c (l;L), on the other hand, does not suffer from 
this complication. By construction at least one site per patch is active, (p) c (/; L) > l~ d , 
so that limx^oo (p) c (I; L) > l~ d . The latter limit is the thermodynamic limit of a 
density and its existence is the most basic assumption in statistical mechanics. Provided 
I is sufficiently large compared to the fixed lower cutoff Iq, so that (J2J) applies, the 
limit implies that lim x ^QQ c (x) > 0, i.e. Q c (x) converges to a non-zero value. This is 
numerically confirmed by the data collapse of (p) c (I; L) in Figured! 

3.2. Variance of the order parameter 

The considerations are similar for the second cumulant, cr^ c (p) (l;L). At I = 1 the 
second moment coincides with the first, so that cr^ c (p) (1 = 1;L) = (p) uc (l — 1,L) — 
(p)uc (I = Ij-^)- I* 1 the unconditional ensemble, this quantity scales asymptotically like 
jj-P/vi. because it is dominated by (p) u (I = 1,L) = (p) u (l,L) oc L^ f3 ^ u± . On the other 
hand, for I = L standard finite size scaling is to be recovered, o\ (p) (/ = L; L) oc L~ 2f3 ^ u± , 
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Figure 2. Data collapse of the variance of the order parameter in the unconditional 
(open symbols) and the conditional ensemble (filled symbols) , as in Figure [TJ The 
data is shown in rescaled form, as a 2 (p) (I; L)l 2 P/ u± vs. l/L. Again, the horizontal 
dashed line is the likely asymptote of the scaling function T c (see ((3])) of the conditional 
ensemble, while the sloped dashed line has slope (3/v±_ ~ 0.795, expected to be the 
asymptote of the unconditional ensemble. The dotted line (apparent asymptote) with 
slope 0.681, however, fits the numerical data very convincingly, indicating that the 
asymptotic regime has not yet been reached. 

so that the scaling function J- U {x) in ([3]) must somehow join these two scaling regimes, 
which implies J-'u(x) = J r (l)x /3 ^ u± for small arguments x. 

However, strictly this argument does not apply, because I = 1 cannot be expected 
to be large compared to the lower cutoff Iq (and in fact is not in the systems 
studied numerically in this article). On the other hand, one might argue that 
a\ (p) (l;L)/a^(p) (/ = 1; L) can be expected to remain finite in the thermodynamic 
limit. While this is not a necessity, the alternative would imply a rather exotic 
behaviour of the variance, with the dotted line (the "apparent asymptote" ) in Figure [2] 
moving further and further away from o~\ (p) (I = 1;L) (which necessarily scales like 
L~P/ V± , dashed line in Figure EJ) with increasing L. If o\ (p) (I; L)/a^ (p) (I = 1;L) 
converges and does not asymptotically vanish in L, then a\ (p) (I; L) inherits the scaling 
of a 2 u (p) (I = 1;L) in L and JF u (x) = F u {l)xPl v± for small x, so that o 2 u (p) (I; L) = 

D u F u ( 1 )l~' 3/v±L ~' 3/ " ± i for sma11 l / L - 

Because the scaling function is a power law only in the asymptote, for intermediate 

values of x = l/L the apparent scaling of a^(p)(l;L) might produce very different 
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effective exponents. This can be seen in Figure H] where the slope suggests J- U {x) oc x 681 . 
A direct estimate of the scaling of a\ (p) (/; L) in Z, at a given, fixed L would then suggest 
a u(p) {h-L) oc W"x +0.681) _ T^jg j s a reminder that scaling assumptions like ([3]) can 
numerically be verified only by a data collapse. The effective exponent of 0.681 is 
of course not a universal quantity and its deviation from 0.795 simply indicates that 
asymptotia has not been reached. 

Not much can be said about the variance in the conditional ensemble. While the 
second moment has a lower bound (namely (p 2 ) c > l~ d ), indicating that its scaling 
function does not vanish in the limit of small arguments, no lower bound exists for 
the variance; in fact a 2 (p) (I — 1, L) — by construction. It is, however, reasonable 
to assume that the variance a 2 (p) (I, L) is finite in the thermodynamic limit for fixed 
/, because by construction every patch always retains some activity regardless of the 
system size. In contrast, in the unconditional ensemble, the moments of the activity 
within a finite fraction of patches might vanish for a duration which increases with 
increasing system size, because for samples to continuously contribute to the average 
only one site needs to be active somewhere in the system. The scaling function T c (x) 
for the conditional ensemble being asymptotically finite, lim^o J~c(%) > 0, is in line 
with the numerical evidence, see Figure [2j 

3.3. Active fraction a(l; L) 

The scaling of the various observables is linked by (00). The fraction of active blocks, 
a(l; L), is given by any ratio of moments taken in the unconditional and the conditional 
ensemble, so that based on the first moments, its scaling is given by 



see (T5]) and (j3J). As can be seen from Equation ([I]), if a moment in the 
conditional ensemble scales like l~ nf3 ^ u , in the unconditional ensemble it will scale like 
l-{n-i)P/u L -p/v_ Similarly, if cr 2 (p) = DJF{l)l~ p l v 'L~^ u for small l/L, the variance in 
the conditional ensemble is dominated by a term proportional to (//L) _2/3//l/ . 

One could consider the fraction of active blocks a(Z; L) as the coarse-grained order 
parameter within a real-space renormalisation group scheme [20], so that (L/l) d is the 
number of coarse-grained sites. The scaling of such an order parameter is proportional 
to (L/l)-Pl u , consistent with 

3.4- Moment ratio R 

The numerical results shown in Figure [1] and Figure [2] all are for the two-dimensional 
contact process. Based on FSS in equilibrium critical phenomena, one would 
normally expect appropriate moment ratios of the (absolute) order parameter, such 
as {p) 2 uc (I; L)/ (p 2 ) uc (I; L) to be universal functions of l/L for I ^> l and to converge 




for / > Iq 
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to a finite value for l/L — > 0. Using Eqs. ([2]) and ([3]), the moment ratio 

^ c(z ' l) = <p 2 u^) (6) 

is universal assuming universality of the scaling functions and of the amplitude ratios. 
However, from what has been said earlier, in the unconditional ensemble, R u (l; L) oc 
(l/L)P^ u± while in the conditional ensemble R c (l] L) indeed converges to a finite value. 

Because R c (l] L) is universal, different models belonging to the same universality 
class, such as the CP and DP, should produce the same values. That is indeed the 
case, as shown Figure [3j For I = L the moment ratios based on the unconditional and 
conditional ensembles coincide by construction and can be compared to the value of 
0.7543(4) [21] (Table VIII, using R C (L; L) = 1/(1 + K 2 /mj)) found by Dickman and 
Kamphorst Leal da Silva (see the dotted line in Figure [3]). The deviation of R c (l] L) from 
R C {L; L) for decreasing l/L does not mean that the former is not universal, just like one 
would generally expect that the value of the latter depends on various geometrical and 
topoligical properties of the lattice, such as its aspect ratio and the type of boundary 
condition [T2] . 

The numerical situation for R u (l; L) (not shown) remained somewhat unclear. Only 
very few points for l/L close to 1 seem to overlap within the error. The system sizes 
simulated did not allow a firm statement as to whether R u (l; L) is actually universal. 

3.5. Surrogate data 

The proposed method can be put to test and compared to others using surrogate 
data, that is data generated in a computer implementation of, for example, the contact 
process, mimicking real-world data as they would be obtained in physical or biological 
systems. In contrast to the simulation data used above, such data consists of a single 
realisation, as if one was to analyse a satelitte image or field data. 

Figure H] shows the conditional box scaling of the activity (p) c (/; L) for a single 
realisation of a two-dimensional system of linear size L = 128 as a function of the box 
size I. The slope of (p) c (I; L) is compared to l~P/ Ul - (thick dashed line) and fits very 
well. No errorbar can be given for / = L, as there is only one such sample, the remaining 
errorbars are estimated from the variance of the corresponding sample of size (L//) 2 , 
assuming independence^ 

For comparison, Figure [5] is a double logarithmic plot of the two point correlation 
function 

2 

(7) 



C(r)=L- 2 Y^I(r')I(r' + r)- ( 




as a function of the absolute distance |r|. Here J(r) indicates occupation of site r and 
the sums run over all L 2 sites r'. It is taken from the same sample as Figure HJ using 

| Obviously, this amounts to an overestimation of the number of independent samples and clashes with 
the assumption of correlations that give rise the scaling in the first place. 
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Figure 3. The moment ratio R c {l\ L) — (p c ) 2 {l\L)/ (p^,) (/; L) for the two-dimensional 
CP and DP. For I > lo (here lo = 64 was needed for satisfactory results) the values 
coincide for both models, as expected by universality. The dashed line shows the 
expected assymptote for small l/L, whereas l/L = 1 is consistent with the finite size 
scaling result R C (L;L) = 0.7543(4) (dotted line) reported in [21]. 



the translational invariance and the eight-fould symmetry of a square. Nevertheless, 
the data is comparatively noisy. Naive scaling arguments along the lines of equilibrium 
phase transitions [22] suggest an asymptote C(r) oc |r|~ d+2 ~' ? = |r| _2,9 ^ ± . In Figure [5] 
it is shown as a thick dotted line, but this asymptote is compatible with the data only 
in a narrow, noisy intermediate regime. 

The correct scaling behaviour of the correlation function however is C(r) oc | r | /^/^-l ; 
which is found by imposing that C(r)/ (p) u does not scale in L [23] (as it would, for 
example, in equlibrium phase transitions, consistent with [21]). Though still noisy, the 
data for small r shown in Figure [U indeed is much better compatible with a scaling 
exponent of —f3/v± = —0.795(6). 

Box scaling of the order parameter is due to the correlations captured in the two- 
point correlation function, as can be seen directly by deriving the variance a\ (p) (I; L) 
from C(r), 

a 2 u (p) (I; L) = r 2d I d d r ! dV C(r - r') oc (8) 

Jld J Id 

consistent with cr 2 (p) (/; L) oc l~P/ u± as observed earlier. Box scaling therefore can be 




Figure 4. Conditional box scaling of the activity as measured in a single instance 
of the two-dimensional CP as surrogate data. The dashed line indicates the expected 
scaling, which fits the data very well. 



regarded as an elegant form of extracting the correlations from the correlation function. 
From that point of view, the advantage of box scaling over a direct investigation of 
the two-point correlation function is merely down to its simplicity: Box scaling is well 
understood theoretically and easy to implement in an experiment, in field work or in a 
computer simulation. 

As a final remark, the quality of the surrogate data shows a strong time dependence. 
Starting from a randomly occupied lattice, chosing too short an equilibration time leads 
to a lack of correlations, with the system still being dominated by the independent, 
random initialisation. Waiting too long, on the other hand, means that the system is 
likely to be low in activity and just about to die out completely, producing sparse and 
biased results. 

4. Discussion 

A block scaling analysis could provide a practical method to overcome the problem 
of limited availability of data in natural systems and allow the analysis of a natural 
system's scaling behaviour without the need of an entire time series. Block scaling 
effectively is a form of sub-sampling and the analysis utilises the universal scaling with 
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Figure 5. The two point correlation function (see Equation ([7])) as measured in a single 
instance of the two-dimensional CP as surrogate data. The error is underestimated by 
assumed independence of the samples. The dotted lines indicate the expected scaling. 
The steeper one corresponds to what is found in equilibrium critical phenomena, which 
is superseeded by the shallower one. This line, with exponent — 0/v.l captures the 
behaviour for not too large distances quite well. 

changing block size, which is characterised in the present work. As it turns out, in order 
to obtain block scaling as known from equilibrium critical phenomena, one has to use a 
conditional ensemble, where moments of the activity in a block enter the average only 
conditional to non- vanishing activity. The situation in absorbing state phase transitions 
therefore is very different from what is expected from equilibrium critical phenomena, 
where no additional condition is needed in order to obtain box scaling corresponding to 
finite size scaling. 

Box scaling effectively measures the correlations between finite boxes throughout 
the system. In (near-)equilibrium systems, i.e. systems with a Hamiltonian, at the 
critical point, the probability density function of the order parameter is symmetric 
around due to the symmetry of the Hamiltonian. As contributions with opposite 
signs cancel, fluctuations decrease with increasing box size and consequently, the box 
averaged (absolute) order parameter and its variance decline. 

In absorbing state systems, this mechanism does not exist: The local order 
parameter is non-negative and therefore cannot cancel out. In an unconditional ensemble 
the order parameter does not change with block size and fluctuations about the mean 
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are not symmetric. For the order parameter to vanish, its fluctuations must vanish 
as well. Nevertheless, introducing the conditional ensemble restores the behaviour of 
(near-)equilibrium systems. 

Within the conditional ensemble, the scaling functions of the observables considered 
converge to a finite value for small arguments, so that the cumulants scale in I with 
the expected exponents, (p) c oc Z~^ x and a 2 (p) oc /~ 2/3 /^J- for sufficiently small 
l/L. Moreover, the moment ratio (p) 2 c / (p 2 ) c is universal and converges to a non-zero, 
universal value for sufficiently small l/L. 

This is not the case in the unconditional ensemble: The scaling functions are 
power-laws themselves and therefore the moments display an unconventional scaling. 
As a consequence, the average activity is constant in /, while its variance scales like 
a 2 (p) oc l~Pl v ^L~~~PI Vl - for small l/L. In addition, the moment ratio (p) 2 u / (p 2 ) u vanishes 
asymptotically in small / / L. The unconventional scaling of the unconditional ensemble is 
due to the existence of the translational symmetry which is not present in the conditional 
ensemble. This symmetry causes the lack of scaling of the first moment, which is 
connected to that of all other moments through (TjQ): 

The ratio of any unconditional moment and its conditional counterpart is the 
fraction of active blocks a(l; L); this quantity itself displays universal behaviour in the 
ratio l/L. Any unconditional moment can be derived from the corresponding conditional 
one and vice versa, by multiplying and dividing by a(Z; L) respectively. 

The key advantage of the conditional ensemble is the scaling of the first moment 
(p) c , which is usually the easiest to determine, carrying the smallest statistical error. 
In the unconditional ensemble, the first moment does not scale at all. Moreover, the 
moment ratio (p) 2 c / (p 2 ) c converges to a finite, universal value, which again is a quantity 
with a comparatively small statistical error. The hope is to use these observables in 
experimental situations where the system is guaranteed to be at the critical point. 
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